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Abstract. In this paper we investigate the impact of delayed tax 
revenues on the fiscal policy out-comes. Choosing the delay as a bi- 
furcation parameter we study the direction and the stability of the 
bifurcating periodic solutions. With respect to the delay we show 
when the system is stable. Some numerical examples are finally given 
for justifying the theoretical results. 
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1. Introduction 



The differential equations with time delay play an important role for econ- 
omy, engineering, biology and social sciences, because a lot of phenomena are 
described with their help. In this paper, we consider a model from economy, 
of the IS-LM type with time delay and we study how the delay affects the 
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macrocconomic stability. The Hopf bifurcation and normal form theories are 
tools for establishing the existence and the stability of the periodic solutions. 
Similar ideas can be found in [5], [9], [10]. 

In [2] Cesare and Sportelli taking into account the papers [6], [7], [8], 
study a dynamic IS-LM model of the following type: 

y(t)=a[J(F(t),r(t))+^-5(F(t)-T(y(0,F(t-r))))-T(y(0,F(t-r)))] 
M{t)^g-T{Y{t),Y{t-T))), 

with Y as income, I as investment, g as government expenditure (constant), 
S as savings, T as tax revenues, r as rate of interest, L as liquidity, M as real 
money supply and a, {3 as positive constants. The time delay r appears in 
function T: 

T{Y{t), Y{t-T))) = d{l - e)Y{t) + deY{t - r), (1) 

where d G (0, 1) is a common average tax rate and e e (0, 1) is the income 
tax share. 

Based on the papers [2], [1], [11], we consider the following IS-LM model: 

y(t)=a[/(F(t),r(t))+^-5(r(t)-r(y(i),F(t-T)))-T(y(t),r(i-T))] 

r{t)=(3[L{y{t)Af))-M{t)] 
kit) = I{Y{t - r),r{t)) - 6K{t) 
M{t)=g-T{Y{t),Y{t-r)), 

(2) 

with 5 > and with the initial conditions: 

Y{t) = ^{t),te [-r,0],r(0) = ri,M(0) = M,,K{0) = K,. 

In the following analysis we will consider the function T given by (1), the 
investment, the saving and the liquidity of the form: 

I{Y{t),r{t)) = aY{t)'^W{t)-''\ a > 0, ai > 0, > 0, 

S{Y{t) T{Y{f),Y{t - r))) = s{Y{t) - T{Y{t), Y{t - r))), s e (0, 1), 

L{r{t)) = mY{t) + , m > 0, 70 > 0, > 0. 

r{t) - r2 

(3) 
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The paper is organized as follows. In section 2 we investigate the local 
stability of the equilibrium point associated to system (2). Choosing the 
delay as a bifurcation parameter some sufficient conditions for the existence 
of Hopf bifurcation are found. In section 3 there is the main aim of the paper, 
namely the direction, the stability and the period of a limit cycle solution. 
Section 4 gives some numerical simulations which show the existence and the 
nature of the periodic solutions. Finally, some conclusions are given. 

2. The qualitative analysis of system (2). 

Using functions (1) and (3), system (2) becomes: 

Y{t)^a[((s-l)(l-e)d-s)Y{t)+de{s-l)Y{t-T)+aY(t)'^^r{ty'+g] 
rit) = l3[mY{t) + ^^^^^ - M{t)] 

k{t) = aY{t - TY'r{ty^ - 5K{t) 
M{t) ^g-d{l- e)Y{t) - deY{t - r), 

(4) 

a,l3 > 0, ai,Q!2 > 0, m > 0, 7o > 0, r2 > 0, 5 > 0, s e (0,1), d e (0,1), 
(0,1). 

System (4) is a system of equations with time delay. The qualitative 
analysis is done using the methods from [3] . 

The equilibrium point of system (4) has the coordinates YqjTq, Kq, Mq, 
where: 

Y,^i [iiizifly-.]-*. A-„ = £(i:_!?2», M.= ™y„+^. (5) 

da rQ — r2 

Using the translation 

F = + Fo, r = X2 + ro, /sT = xa + Kq, M = X4 + Mq 

in (4) and considering the Taylor expansion of the right members from (4) 
until the third order, we have: 

x{t) = Ax{t) + Bx{t - r) + F{x{t),x{t - r)) (6) 
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where x{t) — {xi{t),X2{t),X3{t),X4^{t))'^ , and 
ai = (s — i)d{l — e) — s + apio, bi = {s — l)de, 



apio, hi 



(7) 



-d{l — e), 64 



-de 



ai-l„-a2 
^0 



and 



F{x{t),x{t - r)) = a:(i - r)), F\x{t),x{t - r)), F%x{t),x{t - r)), 

(8) 

F^(x(t),a;(f-r))f, 



F^{x{t),x{t - r)) = —p2oxl{t) + aapiixi{t)x2{t) + —po2xl{t)+ 

YiP^oxlit) + ?>P2lXl{tfx2{t) + ?>Pl2Xi{t)x2{tf + PQzX2{tf) 

F\xit),x{t - r)) = ^i2X2itf + 1730^2(0' 

x(t - r)) = |p2o2;i(i - + apiiXi{t - T)x2{t) + ^PQ2X2{tf+ 

+ '^{pmx\{t - r) + 3p2ia:i(t - Tfx2{t) + 3pi2a;i(t - T)x2{tf + po3a;2(0'^) 

F^(a;(t),a;(t-T)) = 
and 

P20 = - i)i;r^-Vo-"% pii = -«i«2n"'-Vo-("^+'^ 

P02 = «2(a2 + l)yo"^ro-^"^+'\ P30 = «i(q;i - - 2)Y^'-h^^^ 

P21 = -aia2(«i - l)yo"^-Vo P12 = «i«2(«2 + l)yo"^-Vo 

Pos = -a2(«2 + 1)(«2 + 2)yo°Vo-("^+'\ 
12 6 



7i = -7 T^'72 = 7 = 



(ro - r2)2 ' (ro - r2Y ' (ro - r2)^ 



The characteristic equation of hnear part from (6) is: 

det{XI -A- Be~^^) = (A + 5)A{X, r) = (9) 

where 

A(A,r) = P(A) + e-^"g(A) (10) 
P(A) = A=^ + paA^ + piX + po, (5(A) = gaA^ + qiX + go 

and 

P2 = -{aai + /37o7i), Pi = "/?(7o«i7i + mapoi), Po = aapa^poi 
q^^-abi, gi = q;/?7o6i7i, go = a/^Voi- 

To investigate the local stability of the equilibrium point, we begin by 
considering, as usual, the case without delay {r — ). In this case the 
characteristic polynomial is: 

(A + 5)(P(A) + g(A)) = 

hence, according to the Hurwitz criterion, the equilibrium point is stable if 
and only if: 

P2 + '72>0, (pi + gi)(p2 + 52) > Po + go- 
When T > 0, standard results on stability of systems of delay differential 
equations postulate that a equilibrium point is asymptotically stable if an 
only if all roots of equation (10) have a negative real part. It is well known 
that equation (10) is a transcendental equation which has an infinite number 
of complex roots and the some possible roots with positive real part are finite 
in number. 

We want to obtain the values To such that the equilibrium point (5) 
changes from local asymptotic stability to instability or vice versa. We need 
the imaginary solutions of equation A(A, r) = 0. Let A = ±iu; be these 
solutions and without loss of generality we assume c^; > 0. We suppose that 
P{iuj) + Qiiuj) 7^ 0, for all w £ IR. The previous conditions are equivalent to 

(pi + gi)(P2 + g2) T^Po + go- 

A necessary condition to have a; as a solution of A(^a;, r) = is that uo 
must be a root of the following equation: 

/(a;) = uj^ + apio^ + bpco^ + cp = (11) 
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where ap = pi - qi - 2pi, bp = 2gog2 - '2poP2 - qj + pI, cf = pI - go- 
Let k = and = \[f{k)r + l[f'{k)r . 
Using the resuhs from [2], it results: 

Proposition 1. 1. Let ^ > 2' -^^^^ the following cases can be discerned: 
(i) If ap > or bp < 0, then equation (11) has only one real positive 
root; 

(a) If ap < bp > Q and fo < 0, then equation (11) has only one real 
positive root. 

2. If e < ap < 0, bp < and /d < then equation (11) has two real 
positive roots which are distinct if fo 7^ 0. 

3. If e — - and: 

(i) ap < 0, bp — then equation (11) has only one real positive root; 

(ii) ap < 0, bp > 0, then equation (11) has two real positive roots; 

(iii) bp < 0, then equation (11) has only one real positive root. 
Also, we have: 

Theorem 1. If we suppose that the equilibrium point {Yo,ro, Kq, Mq) is 
locally asymptotically stable without time delay, then in conditions of Propo- 
sition 1 there exists only one stability switch. 

Theorem 2. //tq is a stability switch and fo 7^ 0, then a Hopf bifurcation 
occurs at Tq, where 

= —arctg { ^"(^o^^ ~ ~ ^^^^ + ^^^^^ + ^"^^ ~ ^"^i^ 
° Wo V ^0 (91 - 92P2) + i^l{qm - qiPi + pm) - Poqo 

and uo is a root of (11). 

3. The normal form for system (4). Cyclical behavior. 



In this section we describe the direction, stability and the period of the 
bifurcating periodic solutions of system (4). The method we use is based 
on the normal form theory and the center manifold theorem introduced by 
Hassard [4] . Taking into account the previous section, if r = tq then all roots 
of equation (9) other than izcJo have negative real parts, and any root of 
equation (9) of the form A(r) = a(r) ± za;(r) satisfies a{To) = 0, u!{to) = loq 
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and '^'^f^'"''^ 7^ 0. For notational convenience let r — Tq + /i, /J, & H. Then 
dr 

= is the Hopf bifurcation value for equations (4). 
Define the space of continuous real- valued functions as C = C{[—Tq, 0], IR^). 
In T = To + /U, e IR, we regard /j, as the bifurcation parameter. For 
$ e C we define a linear operator: 

L{ii)^ = A$(0) + S$(-t) 

where A and B are given by (7) and a nonhnear operator $) = F($(0), 
r)), where F($(0), r)) is given by (8). By the Riesz representa- 
tion theorem, there exists a matrix whose components are bounded variation 
functions, r]{0,^) with 9 e [— ro,0] such that: 





For $ e Ci([-To,0],IEl^) we define: 



de 



[-To,0) 



/ dr]{t,ii)<P{t), e = o, 

-TO 

0, ^?G[-ro,0) 
F(/.,$), ^? = 0. 
We can rewrite (6) in the following vector form: 

iit = AiiJ,)ut + R{ij)ut (12) 

where u = (mi, M2, u^)'^ , Uf = u(t + 6) for 9 G [— Tq, 0]. 

For ^ e C'^llO, To], IR*^), we define the adjunct operator A* of A by: 



A*'^{s) 



d^iis) 



7 



We define the following bilinear form: 



< ^{9), m >= *^(o)$(o) - / / ^^(e - e)dr]{e)<f>{Od^, 

J -TO J^=0 

where rj{6) = ri{9,0). 

We assume that ±iuo are eigenvalues of ^(0). Thus, they are also eigen- 
values of A*. We can easily obtain: 

^e)^ve^'', ee[-To,0] (13) 

where v = ('i^i, ^2, V3, v^Y , 

(3iai + 64e-^i"'o - mAi) 04 + 646-^1^° 

= 1, V2 = r— T ^ , W4 = 7 , 

Ai(Ai - /37o7i) Ai 

Ai + o Ai(/?7o7i - Ai) 

is the eigenvector of ^(0) corresponding to Ai = ituo and 

*(s)='u;e^2', se[0,oo) 

where w — {wi, W2,W2, W4), 

A2 - Ploli 1 1 . Pl 

Wi = -,W2 = -, W3 = 0,^4 = 

aapoi 1] T] A2?7 

Ai - /37o7i ,1 , , AiToe-^i^o-e-^i^« + 1, /? , . Airoe-^^"" -e"^!^" 

?7= (l+a&i -2 )+V2- — {v4+b4 -2 

aapoi Xf Ai Af 

is the eigenvector of ^(0) corresponding to A2 = —iujQ. 

We can verify that: < *(s),$(s) >= 1, < ^'(s),$(s) >=< ^(s),<l>(s) >= 
0, < ^{s),^s) >= 1. 

Using the approach of Hassard [4], we next compute the coordinates to 
describe the center manifold Qq at fx = 0. Let Ut — Ut{t+ 9), 9 G [— To,0) be 
the solution of equation (12) when // = and 

z{t) =< Ut >, w{t, 9) = ut{9) - 2Re{z{t)^{9)}. 

On the center manifold Qq, we have: 

w{t,9) = w{z{t),z{t),9) 
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where 



w{z, z, 9) = W2o{0)— + wu{e)zz + Wo2{9)— + W3o{9)— + ... 

in which z and z are local coordinates for the center manifold Qq in the 
direction of ^ and ^ and wq2{6) = 'W2o{9). Note that w and Ut are real. 
For solution Ut G flo of equation (12), since = 0, we have: 

z{t)^Xiz{t)+g{z,z) (14) 

where 

g{z, z) = ^>{Q)F{w{z{t), z{t), 0) + 2i?e(z(t)$(0))) = 

= 920^ + ^11^(^)^(^) + 902^ + ^21^^^ + ■ ■ ■ 

where 

520 = ^20^1 + FiQW2, gu = -Fii^i + F^^W2, go2 = F^2'^i + F^2W2, (15) 
with 

= aa(p2o + 2piiV2 + po2V^),Fl^ = aa{p2o + pn{v2 + ^^2) + PQ2V2V2), 

F20 = l0l2vi, F^^ = -fol2V2V2, F^2 = F^o, -^02 = F^, 



and 



where 



g2i = F21W1 + F21W2 (16) 



= aap2o(2^i(0) + ^o(O)) + aapu{2wl^{0) + ?«|o(0) + 2m;i\(0)v2+ 
^0(0)^2) + aapo2{2wf-^{0)v2 + w|o(0)v2) + Q;a(p3o + 2p2i^^2 + '2pi2vl+ 

+P0ZVIV2 + pl2vl + P21^^2) 

-^2^1 = 7o72(2tyii(0)u2 + wio(0)^^2) + 7073^^2^^2- 

The vectors u'2o(^), ■Ji'ii(^) with 9 e [— r, 0] are given by: 

ln\ 920 \,g g02 -X„g , rn 2X^6 

Ai _ 3Ai /"i^yN 

^/3^ Ai6» 511- A26» I r ^ 

■J^ii(o'j = —ve^ — —ve^ + E2 
Ai Ai 
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where 



= -{A + e-^'^°B - 2Ai/)-^F2o, = -{A + 5)-^Fn, 

where F20 = (^20, ^lo, ^lo, Of, Fii = (F^, F^, Ffi , Of . 

Based on the above analysis and calculation, we can see that each gij 
in (15), (16) are determined by the parameters and delay from system (4). 
Thus, we can explicitly compute the following quantities: 

cm = ^(to^n - 2biiP - l\902\') + ^ 

^ "^iA^'^^ ^ '^^ ^ 2Re{Cm)- 

(18) 

and 

y ^ A(g2A^ + giA + go)e-^^ 

3A2 + 2p2\+pi + e-^^(2g2A + gi - T(g2A2 + q^X + qq)) 

In summary, this leads to the following result: 

Theorem 3. In formulas (18), 112 determines the direction of the Hopf 
bifurcation: if 1^2 > 0(< 0), then the Hopf bifurcation is supercritical (sub- 
critical) and the bifurcating periodic solutions exit for r > ro(< Tq); /32 de- 
termines the stability of the bifurcating periodic solutions: the solutions are 
orbitally stable (unstable) if (^2 < 0(> 0); and T2 determines the period of the 
bifurcating periodic solutions: the period increases (decreases) ifT2 > 0(< 0). 



4. Numerical example. 



In this section we find the waveform plots through the formula: 

X(t+^) = ^(t)$(^)+z(t)$(^^^i;2o(^)^'(t)+w^ii(^)-2(t)^(t)H^Wo2(^)^(t)'+Xo, 

where z{t) is the solution of (14), ^{9) is given by (13), W2o{0) , w-ii{9) , Wo2{0) — 
W2q{9) are given by (17) and Xq — {Yq, tq, Kq, Mq)^ is the equilibrium state. 

For the numerical simulations we use Maple 9.5. Wc consider system (4) 
with a = 0.38, a = 0.96, /3 = 1, ai = 0.5, ^2 = 0.83, 70 = 1, d = 0.1, s = 0.3, 
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r2 = 0.003, 6 = 0.2, m = 0.005, g = 50. The equilibrium point is: lo = 500, 
ro = 0.03572181612, Ko = 675, Mo = 33.06065092. In what follows we consider 
two different shares e of delay tax revenues: e = 0.3 and e = 0.8. 

For e = 0.3 we obtain: ^2 = 1.654628706 ■ 10-^ /^a = 2.224294680 ■ 10"^ 
T2 = 2.092652051 ■ 10"^ uq = 0.6685954740, tq = 4.965007916. Then the 
Hopf bifurcation is supercritical, the solutions are orbitally unstable and the 
period of the solution is increasing. The wave plots are given in the following 
figures: 



Fig.l.Waveplot {t,Y{t)) 


Fig.2.Waveplot {t,r{t)) 


Fig.3.Waveplot {t,K{t)) 


mmm- Ij 1 ' 

SOOJOOIfilSS- 

sto- 

4S93999034S- 

msmm- | | | 


1 1 ' M 
1 1 1 U 


3.5721 83 J7E-2 . , , 
3.5721 B31E-2- -, 

3.5721 8236E-2- 
3.5721 81 S3E-2- 
3.5721 8161 E■2■ 
3.67218087E■2■ 
3.5721 SM.? 
3.5721 801 2E-2- U 




6.75000122E2 , ri 
I 

6.75000061E2- 
6.75E2- 

1 

1 

e.74988838E2' i ! 1 , 


1 1 
1 , 




20 40 EO 80 


a SD 40 ED 80 




Fig.4.Waveplot {t,M{t)) 


Fig.5.Waveplot {t,I{t)) 


Fig.e.Waveplot {t,L{t)) 
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Fore = 0.8 we obtain: ^2 = -9.160756314- 10-^ /92 = -1.968119398- 10-^ 
T2 = -1.608068638 - 10-^ ^0 = 0.8553440397, tq = 3.918246696. Then the 
Hopf bifurcation is subcritical, the solutions are orbitally stable and the 
period of the solution is decreasing. The wave plots are given in the following 
figures: 
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Fig.T.Waveplot {t,Y{t)) 


Fig 


.S.Waveplot {t,r{t)) 


Fig.g.Waveplot {t,K{t)) 
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Fig.lO.Waveplot {t,M{t)) 


Fig.ll.Waveplot {t,I{t)) 


Fig.l2.Waveplot {t,L{t)) 
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5. Conclusions. 



From the analysis of the model with continuous time, it results that the 
model accepts a limit cycle. The nature of the limit cycle is given by the 
coefficients (18) which include the parameters of the model. We establish 
the nature of the limit cycle. Because the expressions of the parameter and 
the coefficients (18) are difficult to analyze directly, the use of Maple 9.5 was 
essential. The paper's results confirm that a series of economical processes, 
where a variable with time delay intervenes, have a limit cycle, thus, allowing 
a prediction concerning the evolution of the model. 
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